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We investigate multipartite entanglement in relation to the process of quantum state exchange. 
In particular, we consider such entanglement for a certain pure state involving two groups of N 
trapped atoms. The state, which can be produced via quantum state exchange, is analogous to 
the steady-state intracavity state of the subthreshold optical nondegenerate parametric amplifier. 
We show that, firstly, it possesses some 2iV-way entanglement. Secondly, we place a lower bound 
on the amount of such entanglement in the state using a novel measure called the entanglement of 
minimum bipartite entropy. 



PACS numbers: 03.65.Bz, 03.67.-a 

I. INTRODUCTION 

Multipartite entanglement is entanglement that cru- 
cially involves three or more particles. A well-known ex- 
ample of it occurs in the generalized GHZ state \ip) = 
|0)® M + |1)® M , where M is an integer greater than two. 
Multipartite entanglement is an interesting quantum re- 
source for a number of reasons. Firstly, it is a key re- 
source in quantum computation as i) it has been proved 
that it is a necessary ingredient in order for a quantum 
computation to obtain an exponential speed-up over clas- 
sical computation Q] and ii) it is central to quantum error 
correction which uses it to encode states, to detect errors 
and, ultimately, to help implement fault-tolerant quan- 
tum computation (see, for example, Chapter 10 of 0). 
The second reason why multipartite entanglement is in- 
teresting is that it can manifest nonclassical correlations 
such as GHZ-type correlations 0. Finally, it has been 
conjectured that multipartite entangled states contain a 
wealth of interesting and unexplored physics [j| . 

In order to quantify the amount of multipartite entan- 
glement present in a state, a number of measures have 
been proposed. Firstly, Vedral et al. have suggested a 
measure of multipartite entanglement for a state p which 
is the minimum relative entropy between p and any sepa- 
rable state. For systems with more than two subsystems, 
they defined a separable state as one in which the state 
of at least one subsystem can be factored out from that 
of the others. In addition, Coffman, Kundu and Woot- 
ters |3 have extended the bipartite entanglement mea- 
sure called the tangle to the 3-tangle which measures 
3-way GHZ-type entanglement. Furthermore, Vidal 
has studied entanglement monotones — quantities whose 
magnitudes do not increase on average under local trans- 
formations — and has proposed that all of them can be 
regarded as entanglement measures. Meyer and Wallach 
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[To| have proposed a measure of "global entanglement" 
for n-qubit pure states which is the sum of a number 
of terms involving wedge products. Each wedge prod- 
uct involves computational-basis expansion coefficients 
for various (n — l)-qubit states obtained by deleting a 
qubit from the state of interest. Similarly, Wong and 
Anderson ll] have extended the tangle to an arbitrary 
even number of qubits for pure states. Finally, Biham, 
Nielsen and Osborne |T^ have proposed the Groverian 
entanglement for a pure state based on how success- 
ful Graver's algorithm |13| | performs given the input 
The Groverian entanglement is equivalent to a measure 
in 0, however, it shows an interesting link between an 
entanglement measure and the quantum information pro- 
cessing capability of states. In addition to the measures 
listed above a number of others have also been proposed. 

Quantum state exchange |l4l ITU Il6l Il7l ITsI ] is a newly 
formulated process by which information is transferred 
from an electromagnetic field to the vibrational state of 
one or more trapped atom(s). It is implemented using a 
stimulated Raman process that couples the electromag- 
netic field to the vibrational state and thus transfers in- 
formation from the former to the latter. We explain it in 
more detail in Section ITT1 

In this paper we show that quantum state exchange 
can be used to create an entangled state for 27V trapped 
atoms that is a useful quantum resource. We begin 
in, Subsection III Al by explaining the process of quan- 
tum state exchange via presenting a detailed example 
of it within a simple system consisting of a harmonically 
trapped atom interacting with a cavity mode. Next, Sub- 
section IIIBI shows that quantum state exchange can be 
used to generate an entangled pure state for two groups 
of N trapped atoms located in two spatially separated 
far-off-resonance dipole- force traps (FORTs). In Subsec- 
tion lll Cl we present a detailed summary of the remainder 
of the paper. Sections |lIII and Hvl then investigate the na- 
ture of the state's 27V- way entanglement. That is, the na- 
ture of its entanglement that spans across all 2iV atoms. 
In particular, in Subsections IIII Al and llVl we qualitatively 
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explore this entanglement by presenting a novel neces- 
sary and sufficient condition for the presence of M-way 
entanglement for M-partite pure states and then showing 
that the state satisfies it. In Section Hvl we quantitatively 
explore the state's 27V-way entanglement by first present- 
ing a novel multipartite entanglement measure for pure 
states in Subsection ITV~Al This measure is based on the 
von Neumann entropies |2j, |3j for all the reduced density 
operators obtainable from some pure state of interest by 
tracing over some of the subsystems for the state. After 
defining the measure, we then use it to calculate lower 
bounds on the amount of 27V-way entanglement in the 
state in Subsection IIV Bl Finally, we discuss our results 
m Section El 



II. BACKGROUND THEORY 
A. Quantum state exchange 

Perhaps the simplest system in which quantum state 
exchange can occur [l4| involves a two-level atom con- 
fined within a harmonic trap which, in turn, lies inside 
a linearly damped optical cavity with one lossy mirror 
and one ideal one. The atom's vibrational motion is de- 
scribed by the annihilation operators b x ,b y and b z . The 
two- level atom, which has a transition frequency of Lo a , 
couples to both an intracavity electromagnetic field mode 
of frequency u c described by the annihilation operator a 
and an external laser of frequency lol- The cavity and 
external laser frequencies are chosen so as to drive Ra- 
man transitions that couple adjacent atomic vibrational 
levels. Furthermore, the cavity's axis coincides with the 
x-axis whilst the external laser's beam is perpendicular 
to this axis. Lastly, we assume that the harmonic trap 
is centred on a cavity-field node and thus a schematic 
diagram of this system is as in Fig. 

The system's Hamiltonian is, in a reference frame ro- 
tating at frequency wi, 

H t%t = H sys + + a^R) + H Tes , (1) 

where H rcs is the free Hamiltonian of the reservoir cou- 
pled to the cavity mode, k is a damping rate, R is a 
reservoir operator and H sys is the Hamiltonian for the 
cavity-atom system which is 

H sys = Y 1 j =x ^y yZ hi'j(b^bj + 1/2) + H8a) a + fi,Aer + cr_ 
+h[S L (y, z, t)a+ + £l(y, z, t)aJ] 
+hgo sm(kx) (a V_ + aer+), (2) 

where v x , v y and v z are the harmonic-oscillator frequen- 
cies along the trap's x, y and z axes, cr+ and c_ are 
atomic raising and lowering operators for the two-level 
atom, 5 = lo c — u>l, A = uj a — ujl, El is the complex 
amplitude of the external laser field, k — 2tt/X, where 
A = 2itc/lu c , x — \Jh/2mv x {b x + b x ), where m is the 
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FIG. 1: Schematic diagram of a simple system in which quan- 
tum state exchange can occur. A two-level atom (represented 
by a black circle) lies within a harmonic trap that is itself in- 
side an optical cavity. The cavity, which is aligned along the 
x-axis, supports a mode of frequency ui c and has one lossy 
mirror (with damping constant k) and one ideal one. An 
external laser of frequency ujl is incident from a direction 
perpendicular to the a;-axis. 

mass of the two-level atom, and go (go £ M.) is the cou- 
pling constant for the atom-field interaction. Observe 
that lj c — lol = v x . 

The following reasonable assumptions are made about 
the system so as to make calculations involving it more 
tractable: 

1. The cavity field and external laser frequencies are 
appreciably detuned from tv a and the two-level 
atom is initially in the ground state. Thus, the ex- 
cited internal state is sparsely populated and spon- 
taneous emission effects are negligible and can be 
ignored. 

2. Vibrational decoherence occurs over a timescale 
much longer than that of the interaction producing 
quantum state exchange, as is the case in an ion- 
trap realization of the system |14| . Consequently, it 
has a minimal effect over our timescale of interest 
and is ignored. 

3. The trap dimensions are small compared to the cav- 
ity mode wavelength and thus s'm(kx) <C 1. It fol- 
lows from this that sinfcir ~ Tj x (b x + b x ), where 
r\ x = ky/H/2mi/ x . It also follows that we can ar- 
range things so that the y and z dependence of the 
external laser field is negligible and thus, assuming 
£l is time independent, that z, t) ~ £e~ t ^ L , 
where £ is a real time-independent amplitude. 

4. The damping parameter k is such that i/j>k> 
g r) x £/A. 

Given the assumptions above, and also assuming that 
go /A <C 8 and goia^a) <C £l, H^f° can be rewritten as 
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by adiabatically eliminating the evolution of the internal 
states. Furthermore, for the system under consideration, 
it has been shown |14| that, in the steady-state regime, 
the vibrational state of the atom in the x direction is 
solely determined by the input field (i.e. the light field 
entering the cavity) such that 
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where b(u) = ^ dtb x (t)e^ 
gW x E 2 l{^n) and a in (c) = -fe 

where & in (t) = JJl-oo ^V^""')*^'), where 
Co (a;') is the value of the reservoir annihilation operator 
for the frequency to' at time t = 0. The proportionality 
between b x (uj) and aj n (w) present in Eqn denotes that 
the "statistics of the input field [have been] . . . 'written 
onto' the state of the oscillator" |14| (p. 498) (i.e. onto 
the atom's vibrational state in the x direction). We thus 
say that quantum state exchange has taken place when 
this equation holds. 



B. System of interest 

In this paper we use quantum state exchange to gener- 
ate a particular state involving two groups of N trapped 
atoms which we later show contains multipartite entan- 
glement that is a useful quantum resource. Our work 
follows on from [Tfij in which it was asserted that for 
a certain system consisting of two groups of N trapped 
atoms, quantum state exchange could generate "a highly 
entangled state of all 2[N] atoms". The system we con- 
sider can be seen as a concrete example of that described 
in the last paragraph of [l5j — our main original contri- 
butions are, firstly, demonstrating that the state within 
our system is a multipartite entangled one and, secondly, 
quantitatively analysing the multipartite entanglement 
within it. 

The system we consider comprises of, firstly, a sub- 
threshold nondegenerate optical parametric amplifier 
(NOPA) [H HO, |2l| for which the two external output 
fields first pass through Faraday isolators and then each 
feed into a different linearly damped optical cavity for 
which one mirror is perfect and the other is lossy. The 
axes of both cavities coincide with the x-axis. Each cavity 
supports an electromagnetic field mode of frequency u)j c 
that is described by the annihilation operator a,j, where 
j enumerates the cavities. Within the j th cavity lie N 



identical two-level atoms which each possess an internal 
transition frequency of uj a . These are trapped in a lin- 
ear configuration parallel to the cc-axis by, firstly, a one- 
dimensional far-off-resonance dipole-force trap (FORT) 
[22I |23| . This consists of cavity mode of frequency lot 
which exhibits a standing-wave pattern along the £-axis. 
The frequency lot is strongly detuned from all atomic res- 
onant frequencies and thus the cavity mode's field exerts 
dipole forces on the atoms, trapping each of them near a 
separate node of the field. In addition, the FORT's axis 
is parallel to the x-axis and it thus traps the atoms in the 
x direction. The atoms are also tightly confined in the y 
and z directions by a two-dimensional far-off resonance 
optical lattice [U HiJ [2|| and thus move negligably in 
these directions. The combined effect of all the trapping 
fields is to confine each atom in its own one-dimensional 
trap parallel to the x-axis. Furthermore, the atoms are 
located such that the mean position of each atom coin- 
cides with a node of the cavity field described by a,j . The 

annihilation operator b^ describes the vibrational mo- 
tion of the m th atom in the j th trap in the x direction. 
Finally, external lasers of frequency lul whose beams are 
perpendicular to the x-axis are incident on all atoms and 
thus Fig. illustrates the system under consideration. 
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FIG. 2: Schematic diagram for a system involving, firstly, 
a subthreshold optical nondegenerate parametric amplifier 
(NOPA) whose output modes pass through Faraday isolators 
(represented by an F enclosed in a circle) and then feed into 
linearly damped optical cavities (as indicated by the sinu- 
soidal curves inside both cavites which represent the cavity 
modes ai and 02). These cavities are aligned along the a;-axis 
and both have one ideal mirror and one lossy one (with damp- 
ing constant k). Inside each of them is a far-off-resonance 
dipole-force trap (FORT) that, along with a two-dimensional 
far-off resonance optical lattice, confines iV identical two-level 
atoms in linear chain parallel to the a;-axis. Observe that the 
FORTs' trapping modes are not shown in the diagram. Ex- 
ternal lasers of frequency ujl are incident on each atom in 
both traps from a direction perpendicular to the x-axis. 

The Hamiltonian for the j th optical cavity and the N 
atoms within it is 
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where -ffjo om is the free Hamiltonian for the vibrational 
states of the atoms and Hjq is the free Hamiltonian for 
the cavity field and the atoms' internal states. The term 
Hjj is the interaction Hamiltonian describing the Ra- 
man processes involving the cavity field, the external 
lasers and the atoms. To be more specific, H^° m — 

hvjx El=i t ^™'' + , where v jx is the vibrational 

frequency of the FORT in the j th cavity (which we call 
the j th FORT) in the x direction. This frequency is equal 
to LUj C — lul. The Hamiltonian Hjo is, in a frame rotating 
at frequency ujl, 
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(5) 4. All atoms are tightly confined in the y and z direc- 

tions. This allows us to ignore the y and z depen- 
dence of the external laser fields and thus, assuming 
£l is time independent, it follows that £l{v, z, t) ~ 
g e -i<t>L ( wnere £ is a real time-independent ampli- 
tude. 

5. The damping parameter k is such that Uj X >k> 
goVjx\^N£ I A, where rjj X = kj^Jh/2mvj X , where m 
is the mass of each atom. 

Given the assumptions above, we can write Hj total in 
terms of normal-mode creation and annihilation opera- 
tors (by adiabatically eliminating the evolution of the 
internal states) as 

(6) JTj total = h jr Vjx ( B ]^B^ + I) + hSa^a 



where 5j = uij C — u>l, A = uj a — u>l and crj™ 1 ' and 

are raising and lower operators for the internal states of 

the m th atom in the j th trap. The term Hji is 
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(8) 
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where b[™^ is the annihilation operator for the m th nor- 



hgo sia(kjXj m ){ajOj_ + Ojcr" ), (7) mal mode for the j th trap in the x direction. For exam- 
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where £l is the complex amplitude for all external lasers, 
kj = ujje/c and go (go € R) is the coupling constant for 
the atom-field interaction. Finally, Hj rcs is the Hamil- 
tonian for the external reservoir that couples to the j 
cavity for which Rj is a reservoir annihilation operator 
and k is a damping constant. 

The feasible assumptions below are made about the 
system in order to simplify calculations for it and to focus 
on its most important aspects: 

1. The cavity field and external laser frequencies are 
appreciably detuned from u a and all two-level ions 
are initially in the ground state. Thus, the excited 
internal states are sparsely populated and sponta- 
neous emission effects are negligible and can be ig- 
nored. 

2. Vibrational decoherence occurs over a timescale 
much longer than that of the interactions produc- 
ing quantum state exchange and consequently can 
be ignored. 

3. The wavelength of the cavity mode described by a } - 
is much greater than the distance that any atom in 
the j th trap strays from the cavity-field node about 
which it is trapped. Thus, sin(fcxj m ) ~ kxj m << 1 
and hence all atoms experience a potential that is, 
to a good approximation, harmonic. This justifies 
the form of 7J^ tom . 



pie, B^J is a centre-of-mass mode annihilation operator 



which is B$ = 1/VN (bfj + bfj + ... whilst Bf] 
is the annihilation operator for the breathing mode which 
is B?J = 1/V2 (-&£> + bfj) when N = 2. 



Comparing Eqn (JHJ to Eqn (J3J, we see that B^} in 
Eqn (JSJ) plays an almost identical role to that of b x in 
Eqn ©. Given that s'm(kx) ~ r] x (b x + b x ) in Eqn 0, 
the only difference between the forms in which the two 
operators appear results from a factor of v^V" appear- 
ing in front of £ in Eqn JHJ. As a consequence, B^ in 
Eqn JHJ couples to the cavity mode dj identically — aside 
from the factor of V^/V — to the manner in which b x cou- 
ples to a. It follows that as quantum state exchange takes 
place in the system described by H^^ e with information 
about an input electromagnetic field being transferred to 
b x , it also occurs in the system described by Hj tota i due 
to the correspondence between the two system's Hamil- 
tonians. Thus, in the latter system, information about 
the input field is transferred to the centre-of-mass mode 
for the trapped atoms in the x direction just as if this 
mode was a vibrational mode for a single harmonically 
trapped atom. The only difference between the iV-atom 
case and one described by H^f° is that the effective cou- 
pling in the former case is increased by a factor of y/N. 
This conclusion can also be verified via comparing the 
Langevin equations for and b x . Due to symmetry 

considerations, collective modes other than the centre- 
of-mass modes do not absorb any photons in modes a\ 



5 



and ci2- Thus, assuming these other modes are initially 
in vacuum states, they remain so during quantum state 
exchange. 

In |15|. it was shown that we can transfer the intra- 
cavity steady-state for the subthreshold nondegenerate 
parametric amplifier which is 

|,/;> = coi7 s ^ otanhnr|n>i|n>2 ' (9) 

where the subscripts 1 and 2 denote the two output 
modes and r is a real squeezing parameter, into the vi- 
brational states in the x direction for two single trapped 
atoms in different harmonic traps. Using the correspon- 
dence between the quantum state exchange processes in- 
volving a single harmonically trapped atom and N har- 
monically trapped atoms demonstrated above, it follows 
that in the system illustrated in Fig. [21 we can transfer 
into the centre-of-mass modes in the x direction for 
the two sets of N trapped atoms thus producing, in the 
steady state, 

IV>cm> = ^ 7 S^ =0 tanh AA r|AA) 1 |AA) 2 , (10) 

where \Af)j denotes the centre-of-mass vibrational num- 
ber state for the x direction with eigenvalue AT for the 
atoms in the j FORT. In writing this state we have 
omitted the states of collective modes other than the 
centre-of-mass modes as we have assumed these other 
modes are in vacuum states throughout the quantum 
state exchange process. 

Importantly, the process of creating IV'cm) just out- 
lined does not seem to be overly experimentally infea- 
sible. This is so as optical cavities and nondegenerate 
optical parametric amplifiers have been widely realized 
quantum optical laboratories for some time. In addition, 
neutral atoms have been confined within standing-wave 
dipole- force traps that, in turn, lie within optical cavities 
[29j. Relatedly, experiments in which a single harmoni- 
cally trapped ion has been placed within an optical cavity 
have been conducted [28J. 



C. Summary 

In this paper, we explore multipartite entanglement 
in relation to quantum state exchange and in Section fllll 
follow on from a multipartite entanglement condition im- 
plicit in work by Diir and Cirac |30| by presenting a novel 
condition. The satisfaction of this novel condition implies 
that any pure state comprising of M subsystems is M- 
way entangled. Here, an M-way entangled state is one 
possessing entanglement that spans across M subsystems 
as does the generalized GHZ state |V>) = |0)® M + |1)® M . 
After presenting this condition, we then use it to show 
qualitatively that IV'cm) is 2iV-way entangled. In Sec- 
tion IIVI we quantitatively consider the entanglement in 
IV'cm)- We introduce a novel multipartite entanglement 



measure for pure states we call the entanglement of min- 
imum bipartite entropy or £mbe which is the minimum 
of the von Neumann entropies of all the reduced den- 
sity operators obtainable from some pure state of inter- 
est by tracing over some of the subsystems for the state. 
After this, we use £mbe to calculate a lower bound for 
the amount of 4-way, 6-way and 8-way entanglement in 
IV'cm) for N = 2, 3, 4 respectively for a range of r values. 
Finally, we discuss the nature of our results. 

It is interesting to investigate the nature |V»cm)'s 2N- 
way entanglement for a number of reasons. Firstly, it 
has been claimed — but not demonstrated — IV'cm) is 
"entangled state of all 2[N] ...atoms" ^4|- ft is thus 
interesting to investigate |^cm)'s 2iV-way entanglement 
in order to see if this implied claim is true. Secondly, 
it is interesting to investigate IV'cm) 's 2iV-way entangle- 
ment as it is a massive-particle state which is important 
as, to date, mostly massless photons have been used to 
experimentally investigate entanglement. Thirdly, if the 
claim is true, then it means that IV'cm) is a state consist- 
ing of 2N entangled harmonic oscillators, each possessing 
an infinite- dimensional Hilbert space as opposed to the 
two-dimensional Hilbert space of a qubit, that is 2N-w&y 
entangled. 



III. QUALITATIVE RESULTS 
A. Negative partial transpose sufficient condition 

Assume that, for a certain state p, we wish to know the 
answer to the question "Does p contain at least some M- 
way entanglement?" Whilst answering this question does 
not tell us everything about the nature of p's M -way en- 
tanglement, it nevertheless tells us something of interest. 
One way to answer it, provided that p consists of qubits, 
is to use a condition that can be readily derived from 
work by Diir and Cirac [30| . This condition involves neg- 
ative partial transposes (NPTs) |3lJ, |32|, |33J and thus we 
name it the NPT sufficient condition. It is sufficient for 
the presence of M-way entanglement for all p's consisting 
of P qubits, where P > M, and is based on generalizing 
the notions of separability and inseparability to many- 
qubit systems. Before stating the condition, it is first use- 
ful to mention two things. Firstly, we define an M-partite 
split of p |3(il ] to be a division or split of p into M parts 
which each consist of one or more subsystems. Secondly, 
we observe that p can always be converted to a state that 
is diagonal in a certain basis by a "depolarization" pro- 
cess consisting of particular local operations 30] . This 
basis consists of M-qubit generalized GHZ states of the 
form \tp) = 1/V2(|j)|0) ± |2 JV " ;L - j - 1)|1)), where j is 
a natural number that we write in binary as M — 1 bits, 
i.e. j = jij 2 ■ ■ -Jm-I) where j x is the x th bit in fs bi- 
nary representation. Given these two things, the NPT 
sufficient condition states that p is M-way entangled for 
a given M-partite split if the diagonal state that it de- 
polarizes to is such that all bipartite splits that contain 



6 



the M-partite split have negative partial transposes. A 
bipartite split is one that divides a system into two parts, 
i.e. a 2-partite split. Also, a bipartite split that contains 
an M-partite split is one that does not separate members 
of any of the M subsystems onto two different sides of 
the bipartite split. That is, one that does not cross any 
of the divisions created by the M-partite split. 



B. Result 

Following on from the NPT sufficient condition, we 
propose a necessary and sufficient condition for the exis- 
tence of M-way entanglement for M-partite pure states. 
Our condition is based on the traces of the squares of re- 
duced density matrices obtained by tracing over some of 
the subsystems constituting our system of interest. Af- 
ter formulating it, we then use it to demonstrate that 
IV'cm) contains some 2iV-way entanglement. Our moti- 
vations for employing our condition instead of the NPT 
sufficient condition are that i) it seems to be mathemat- 
ically simpler to calculate whether or not our condition 
is satisfied and ii) as we are concerned with a pure state, 
our condition is stronger than the NPT sufficient condi- 
tion in the sense that it is both necessary and sufficient 
as opposed to just being sufficient. 

Our M-way entanglement condition utilizes the fact 
that when a pure state \ip) for M subsystems is M-way 
entangled then we cannot write it as \ip) = |0i)Q j <8> 
\ < h)Q-> where and \4>2}q- are the states for the 

subsystems denoted by Qj and Qj respectively and both 
Qj and Qj denote at least one subsystem. To put this 
another way, when \tp) is M-way entangled then there is 
no way to represent it as the tensor product of two pure 
states. Consequently, excluding all such possibilities suf- 
fices to show, and is also, in general, necessary to show, 
that |0) is M-way entangled. This can be done by first 
checking that no single-subsystem state can be factored 
out from the state of the remaining M — 1 subsystems. 
We do this by checking that the traces of the squares 
of all the reduced density operators obtainable from \xp) 
by tracing over one subsystem are less than one. That 
is, that TrdpQj} 2 ) < 1, where pQ j is the reduced den- 
sity operator obtained from |0)(0| by tracing over the 
subsystem denoted by Qj for all Qj denoting just one 
subsystem. We can then repeat this procedure, consid- 
ering all Qj's corresponding to all pairs of subsystems, 
then all triples and so forth until we have considered all 
Qj's corresponding to all sets of R subsystems, where 
R = [_M/2\ , where [x\ is the largest integer less than or 
equal to x. It is sufficient to only consider sets of up to 
those corresponding to *lM/2\ subsystems as a necessary 
condition for being able to factor out any larger a number 
of subsystems from \ip) is the ability to also factor out 
[M/2\ or fewer subsystems. Underlying the process just 
described is that of seeing whether or not we can exclude 
all the ways that |0) could fail to be M-way entangled. 

Our condition can be formalized as Definition 1 



which is as follows: 

Definition 1 : For a pure state \tp) for M 
subsystems, consider the set Q whose mem- 
bers Qj are themselves sets of subsystems 
for the system corresponding to \ip). This 
set Q contains all sets of P subsystems for 
this system, where 1 < P < [M/2\. Given 
this, \ip) is M-way entangled iff, for all Qj, 
Tr([ / 9q i ] 2 ) < 1, where pQ j is the reduced 
density operator obtained by beginning with 
an d tracing over the subsystems Qj. 

To illustrate Definition 1 consider, for example, the 
GHZ state |V>)ghz - l/\/2 (|000)i 23 + |Hl)i2 3 ), where 
the subscripts 1, 2 and 3 denote subsystems of |0)ghz- 
The parameter P = [3/2\ = 1 and consequently the 
set Q comprises of all sets of one subsystem and thus 
Q = {{1}, {2}, {3}}, where the numbers again denote 
subsystems for |0)ghz- For the element {1}, for exam- 
ple, Tr([p {1} } 2 ) = 1/2. Calculating Tr([p Q] } 2 ) for all of 
Q's other elements, we find that it is 1/2 in all three 
cases. Thus, |0)ghz satisfies Definition 1 and hence is 
said to be 3-way entangled, as is the case. 

To further explain Definition 1, we now apply it to 
determining whether the following four-party states are 
4-way entangled : 

1. ) - l/\/2(|0000)i2 34 + |1111) 1234), 

2. ) |^ 2) ) = l/v^|0)i ® (|000) 234 + |1H)234) 

3. ) |4 3) ) = |0+) 12 ® |0 + ) 34 . 

Turning to 1.), we see that upon tracing over any sin- 
gle subsystem, we produce a reduced density operator of 
the form p Qj = 1/2(|000)(000| + |111)(111|) for which 
TVQpqJ 2 ) = 1/2. Similarly, tracing over any two sub- 
systems produces a density operator of the form pq . = 
1/2(|00)(00| + |11)(11|) for which, again, Tr{[p Qj ] 2 ) = 
1/2. Thus, Definition 1 gives the correct result that 

|04^) is 4-way entangled. For 2.), tracing over the first 
subsystem produces |^ 3) ) = 1/^(1000)234 + |111) 2 3 4 ), 
which is a pure state and hence Tr^pQ^ 2 ) = 1 for the 
corresponding j. Consequently, Definition 1 tells us 

that |V ; 4 2 ' ) ) is not 4-way entangled, as is the case. For 3.), 
tracing over any one subsystem produces the mixed state 
p = I/2(^\<j) + } 3 4{(j) + \ and so we might be tempted to infer 

that I ■04 ) is 4-way entangled. However, when we trace 
over subsystems 1 and 2 or subsystems 3 and 4 we pro- 
duce the pure state |0+) for which Tr ([|0+)(0+|] 2 ) = 1. 

Hence, Definition 1 correctly tells us that \ip^) is not 
4-way entangled. 

In applying Definition 1 to |0cm), we first write 
I 0cm) in terms of vibrational number states for the 2N 
atoms involved as we wish to see if they are 2A^-way en- 
tangled. As a step towards doing so, upon observing 

th- IA0, = lM |0),. we express Wl in 
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terms of vibrational number states in the x direction for 
individual atoms as 



a.c.(n,JV") 



(ii) 



a 



where n is the iV-component vector (ni,ri2 ■ ■ - tin), the 
state \ft)j = \ni)j ® |t&2)j ■ ■ • |^jv)j) where |rife)j denotes 
number state for the k th atom in the j th FORT and 

c(n,M) = j{n\N)j 

r N\fN-nx\ M- rii... - njy-2 



VW! x 



x y 



ni!ri2!...njv! 



(12) 



The sum c ^ denotes the sum over combina- 
tions of ni, «2 . . . n^v such that X^jLi n j = N [35j]. Using 
Eqn 1111) to represent IV'Cm) in terms of vibrational num- 
ber states for individual atoms, we obtain 



1 



cm; 



coshr 



tanh 



N=0 



r i E c(U,J\f)\n) 

a.c.(n,Af) 



E c{fh,M)\fa)2 

, a.c.(m,AA) 



(13) 



We now show that the right-hand side of Eqn l|13[) sat- 
isfies Definition 1 and thus that I^Cm) is 2iV-way entan- 
gled. We do this by first writing \ipcM) as the most gen- 
eral bipartite state possible involving vibrational number 
states for individual atoms. Next, we show that, upon 
tracing over the atoms in the half of the bipartite split 
containing the lesser number of atoms and then finding 
the trace of the square of the resulting reduced density 
operator, that this is less than one. It follows that, for all 
j, 7 1 r([pQ i ] 2 ) < 1. Hence, we satisfy Definition 1 and 
so I^Cm) is 2A-way entangled. 

Dividing the atoms in IV'cm) into subsystems A and B 
containing, respectively, R and 2 N — R atoms (R ^ 0), 
we can write I^Cm) as 



oo 

\ipcu) = y ^2 l c i \fi) A 

t=0 



\9t)B, 



(14) 



where = HIsObII = 1 and the \fi) A , but not nec- 

essarily the \gi)Bi are mutually orthogonal. (As we can 
always write \iPcm) in biorthogonal form j3^|, there exist 
\gi)B that are mutually orthogonal. However, we are not 
concerned with this form in the current calculation and 
so do not consider such a decomposition of IV'Cm)-) To 
give an example, when N — 2 and A contains the first 
atom in the first trap 



I^cm) 



1 



coshr 



|o>a® |ooo) 



tanh ) 



|101)i 



tanhr 



110), 



tanh 2 r 



|202), 



v^tanh 2 



1211) 



tanh r 



|220), 



1 , , / tanh r , 

Da® ( — ^— |001} B 



coshr 
tanh r 



V 2 



|010)j 



tanh" r 



V8 



|102)j 



V2tanh 2 r 



111) 



tanh 2 ; 



1120), 



1 , . / tanh 2 r , 

|2) A ® — -. — |002) B 



coshr 



V2 tanh 2 



V 4 

1011)* + 



tanh r 



|020) B + 



(15) 



where \x)a = \ni = x)a and ^1X2x^)3 = \n>2 — Xi,mi = 
X2,ni2 — xz)b- Here, for example, cq = 1/ coshr, 
\/o)a = |0)a, ci = 1/ coshr, \fi) A = \1) A , 



1 



\9o)b 



/Mo 
tanh r 



1 000) j 



tanh r 



101) 



H0)b|) 



tanh 2 : 



1202), 



v^tanh 2 



1211), 



and 



\9i)b = 



1 / tanh? 



/Mi V 2 



1001), 



tanh 2 ; 
tanhr 



|220) B 4 
|010) fl 



+- 



tanh 



tanh 2 r 



-|102)j 



V^tanh 2 



1111) 



|12Q)j 



where Mo and A^i normalize \go)B and \gi)B- Upon 
tracing over A in Eqn l|14|l and squaring the resulting 
reduced density operator pq a , we obtain 



[PQ A ? = E c i c j\9i){9i\9j){9j\ 

Calculating the trace of [pq a ] 2 yields 

00,00 

Tr{\p QA f)= E c " c H d «| 2 ' 
»,i=o 



(16) 



(17) 



where dij = (gi\gj)- As the trace of a density operator is 
always one, we know that 



E 

i,j=0 



2 2 
CiCj 




(18) 
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It thus follows from Eqn (|17fl that, as Cj 7^ for all i, if 
\dij\ 2 < 1 for at least one dij then Tr ([pQ A ] 2 ) < 1- 

As the centre-of-mass state |7V)i|7V)2 has an even num- 
ber of centre-of-mass phonons in total (27V), when we 
express it as a sum of vibrational number states for in- 
dividual atoms, these states all contain an even num- 
ber of individual phonons in total. Furthermore, because 
\iPcm) contains the state \Af = 0)i|7V = 0)2, one \fi)A in 
Eqn (|T3*jl . which we denote by \ff eTO )A, is a tensor prod- 
uct of ground states for some of the 27V atoms in IV'cm)- 
For example, in Eqn (JI3J, |/f ero )^ = 10)^. Given that, in 
general, |/j Zero )A contains zero individual phonons, only 
states with an even number of individual phonons in to- 
tal are present in the |<?i)s with the same index i, which 
we denote by |gf ero ) b- This is so as we require the total 
number of individual phonons in \ ff oro ) a ® |<?i° ro ) b to be 
even. 

In addition to |/f ero ) J 4, because (i/'Cm) includes the 
term \Af = l)i|A/" = 1)2, there also exists an \fi)A in 
Eqn l|13l) containing just one individual phonon which we 
denote as \f1 nc ) A - For example, in Eqn JTSJ \f? ae ) A = 
\1) A - In general, the \gi)B with the same index i as 
\f° nc )A, which we denote by |g° ne )s, comprises of states 
with an odd number of individual phonons in total as 
dictated by the requirement that the total number of in- 
dividual phonons for |/° i1c )a (8 \g° nc )B is even. Thus, 
|<7° ne )B is orthogonal to |gf ero )B and the corresponding 
|4,.|2 = |^zero| 5 one^|2 = q R eturn i n g to the right-hand 

side of Eqn (fP7|l . this means that Tr ([pq a ] 2 ) < 1 for Qa 
and thus that Definition 1 is satisfied. This allows us to 
infer that IV'Cm) is 27V- way entangled and consequently 
we have verified the assertion that |?/>cm) is an "entan- 
gled state of all 2[N] ... atoms" — except, of course, when 
r = 0. 



IV. QUANTIFYING THE AMOUNT OF 
2 TV- WAY ENTANGLEMENT IN |i/>cm) 

A. Theory 

In the previous subsection we presented a qualitative 
result which showed that IV'cm) possessed some 27V- way 
entanglement. However, we would also like to know how 
much 27V-way entanglement IV'cm) contains. For this 
reason, we present an novel quantitative measure of M- 
way entanglement for M-partite pure states, for arbitrary 
M. This measure is based on the von Neumann entropies 
of reduced density operators produced by considering all 
bipartite splits for some state of interest. We call it 
the entanglement of minimum bipartite entropy or TSmbe, 
which we soon define. After this, we then argue that it is 
a plausible measure and finally use it to calculate a lower 
bound on the amount of 27V-way entanglement in IV'cm)- 

For a pure state with M subsystems, -Embe is 

£mbe(|V>» =mm(S aU ), (19) 
where S a n is the set containing the von Neumann en- 



tropies of all the reduced density operators obtained from 
|'0)(V'l by tracing over a set of P subsystems in where 
1 < P < [M/2\. The function min(A) returns the 
smallest element of the set X. Thus, as the von Neu- 
mann entropies of both sides of any bipartite split of 
\%l>) are equal 0] (p. 513), S a n contains the von Neu- 
mann entropies for all the reduced states that we can 
generate from \ip). For example, when \tp) = |^)ghz = 
1/V2 (|000) 123 + |lll)i23), the sets of subsystems con- 
taining P members that we trace over in obtaining S a n 
are {1}, {2} and {3}, where the numbers denote either 
the '1', '2' or '3' subsystems of |V")ghz- As the von Neu- 
mann entropy of the state p is S(p) = — Tr(p\og 2 p) 0,0, 
the von Neumann entropy for the reduced density op- 
erator generated from |V')ghz('0I upon tracing over the 
subsystem denoted by any one of these sets is 1. Hence 
Sail = {1, 1, 1} and so -EmbeGV^ghz) = 1- We thus we 
say that \4>)gwz has 1 unit of 3- way entanglement. 

To provide some insight into £mbe, it is now shown 
that it can be thought of as a distance-based measure 
of M-way entanglement. That is, as measuring the dis- 
tance between |?/>) and the closest pure state with zero 
M-way entanglement given a certain metric. To under- 
stand this, observe that, naively, it seems reasonable to 
think that there exists a pure state IVWo) with zero M- 
way entanglement that has an identical S a n to |'0)'s ex- 
cept for one element. This element corresponds to the 
smallest element of S a n(|-0)) and is zero. The next step 
in comprehending the distance-based nature of £mbe is 
representing Saiid^)) and S a n(|VWo)) by points A and 
B respectively in a co-ordinate space for which each co- 
ordinate denotes the possible values of an element of ei- 
ther S a n(|'0)) or S a n(|-0 ZG ro))- That is, a space that graph- 
ically represents S a n(|i/')) and S a n(|?/; zcro )). For such a 
space, we observe that no pure state with zero M-way 
entanglement is represented by a point closer to A than 
73. It is in this sense that we think of |f/ , zero) as being the 
closest pure state to \ifj) with zero M-way entanglement. 
Finally, the distance-based nature of 75mbe(|V')) can be 
seen by observing that the distance between A and 73 is 
^mbe(IV'))- This point is illustrated in Fig.|3|for the 3- 
way entangled state \4>) comprising of three subsystems 
for which S a ii(|</>)) = {Si, S 2 , S 3 }, where Si < S 2 ,S 3 
and Si, 52, S3 ^ 0. Naively, the closest pure state to \(f>) 
with no 3- way entanglement \4> ze ro) seems to be such that 
S a u(|<2Wo)) = {0,S 2 ,S 3 }. Representing S a ii(|</>)) and 
S a ii(|</>zero)) graphically in the manner described above 
by points A and B in Fig. |3 we observe that the dis- 
tance between them is Si = 75mbe- Generalizing this 
notion, we see that 75mbe can be viewed as measuring 
the distance between \ip) and the nearest pure state with 
zero M-way entanglement. This distance seems to be a 
plausible measure of |-i/>) 's M-way entanglement and thus 
7?mbe appears to have an underlying intuitive motiva- 
tion. 

To further highlight the plausibility of 75"mbe , consider 
the following analogy. Imagine an ordinary chain with 
M links. If M — 1 of these are strong and the other 
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S 2 

FIG. 3: Co-ordinate space illustrating that -Embe can be 
seen as a distance-based entanglement measure. The distance 
between A (the point representing S a n(|0}) = {Si, 5*2, S3}, 
where \<j>) is 3- way entangled and Si < S2,Ss) and B (the 
point representing S a n(|<?f>zcro)) = {0, S2,S3}, where |0 ze ro) 
appears to be the closest pure state with no 3-way entan- 
glement to \(f>)) is -Embe. The quantities Si, S2 and S3 are 
dimensionless. 

one is weak, then the chain is close to breaking and so 
only has a small amount of "nonbrokcn-ness" — even 
though all but one of the links are solid. This is so as 
nonbroken-ness is a wholistic property that is a manifes- 
tation of the nature of all M links. Relating this to Embe, 
just as nonbroken-ness is a wholistic property, so Embe 
measures a wholistic property, namely M-way entangle- 
ment, that relates to the nature of all M subsystems of 
M-partite states. In analogy with a chain with just one 
weak link, an M -partite pure state for which all members 
of S a ii are large, except for one, is very close to possessing 
no M-way entanglement. In this way, we see that -Embe 
and, in particular, the presence of the min function in it 
seem plausible. 

Another positive feature of -Embe is that it satisfies 
three well-known desiderata for bipartite entanglement 
measures |36j, as we now show. (It seems plausible that 
these should also be desiderata for multipartite entangle- 
ment measures.) They are : 

1. ) The proposed entanglement measure is zero for all 

product states. 

2. ) The proposed entanglement measure is invariant 

under local unitaries. 

3. ) The proposed entanglement measure does not 

increase on average under local operations, clas- 
sical communication (LOCC) and division into 
subensembles. 

Beginning with 1.), if the state of interest is a product 
state, where we define a product state to be one for which 
we can factor out the state of at least one of the subsys- 
tems, then at least one member of S a n is zero and so 



Embe is also zero, as we desire. Turning to 2.), we note 
that for a general bipartite split, the von Neumann en- 
tropy of the reduced density matrix obtained by tracing 
over the subsystems on the side of the split with the lesser 
number of particles is invariant under unitary transfor- 
mations which act on only one subsystem. Consequently, 
if we define local unitaries to be those which act just on 
a single subsystem, then Embe satisfies 2.) 

In considering 3.), it is important to remember that 
Embe is only for pure states and thus we ignore local op- 
erations that convert \ip) to a mixed state. For example, 
we do not consider local operations that transform \ip) 
to a state that is close to a maximally mixed state and 
thus has large values for the von Neumann entropies of 
all its reduced states. We choose this example as such 
local operations increase the value of min(S a n) for a sys- 
tem of interest. However, they manifestly do not increase 
its M-way entanglement but instead transform its state 
into one for which Embe is not applicable. With this 
constraint in mind, we define a local operation to be one 
that involves just one subsystem, such as a projective 
measurement on a single subsystem. Given this defini- 
tion, it can be shown that for bipartite pure states LOCC 
and division into subensembles cannot increase the av- 
erage entanglement of any state as measured by the von 
Neumann entropy of its reduced states (entropy of entan- 
glement) [3(| . It follows that they also cannot increase 
any member of S a n(|?/>)), on average, as these faithfully 
measure the bipartite entanglement in \ tp) given some bi- 
partite split for it. Thus, Embe also cannot increase, 
on average, under LOCC and division into subensembles 
and so Embe satisfies 3.) 

Another well-known desideratum for a bipartite entan- 
glement measure is that it is additive over tensor products 
36] . However, it can be shown that Embe is superaddi- 
tive. That is, that the M-way entanglement of a com- 
bined state generated from two states with a and b units 
of M-way entanglement can be greater than a + b (but, 
importantly, not when M — 2). It is an open question as 
to whether or not multipartite entanglement is additive 
and so we do not know if the superadditivity of Embe 
represents a flaw. 

For Embe to be a reasonable measure, it ought to re- 
duce to the standard pure state bipartite entanglement 
measure of the entropy of entanglement. For Embe, 
when N = 1, we have Embe = niin(S a n) = S{i}, 
where S^ij is the von Neumann entropy for the reduced 
density operator Pq {1} = Tri(|-0) ('(/'I), and so we re- 
cover the desired measure, namely the entropy of en- 
tanglement. Finally, Embe seems to be plausible as for 
= ,Jc\0)® N + yfi~~c\l}® N , where c € [0, 1] and N is 
a positive integer, Embe = — clog 2 c— (1 — c) log 2 (l — c). 
This expression increases monotonically in the interval 
c G [0,1/2] and attains its maximum value of one for 
c = 1/2. Such behaviour seems reasonable. 
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B. Results 

In this subsection we use -Embe to calculate lower 
bounds on the amount of 2iV-way entanglement present 
in \ipcu) f° r N — 2, 3, 4, for a range of r values. We ob- 
tain these lower bounds by, first, calculating TVQ/uqJ 2 ) 
for a general Qj. Next, we determine the linear entropy 
Sl(pq ) Ei2| from the relation S L (p Qj ) = 1 - Tr([p Qj } 2 ) 
and then use the fact that Sl (p) / log 2 e < S(p) to obtain 
our lower bounds. We calculate a lower bound rather 
than -Embe itself as it is computationally infeasible to 
calculate Embe due to the fact that it is computation- 
ally infeasible to calculate the required von Neumann 
entropies of reduced density operators given the infinite- 
dimensional bases of the harmonic oscillators comprising 
\ipCM)- This is so as these are generally calculated by 
first diagonalizing p and it is computationally infeasible 
to do this, in general, when p is a square matrix of infinite 
dimension. 

We begin with the initial density operator pcm = 
IV'cmKV'cmI which can be written in the centre-of-mass 
number-state basis as 

oo,oo 

Pcm - £ /(AA,AA')|A0i|A0 2 2(M'\ 1 (JV'\, (20) 

where f(N,N') — tanh JV ' + r/ cosh 2 r. To obtain a gen- 
eral pq . , we trace over the first T atoms in the first FORT 
and the first V in the second one, arriving at 

66 oo,oo 

PQi =EE f{MM')(PW)iW) 2 2 (AA'| 1 (AA'|P), 

J21) 

where P is a dummy variable given by P = 

(Pi\p2 ] ■■■Pt^Pi' '•••Py )i where pa denotes a vibra- 
tional number state for the a th atom in the x direc- 
tion in the j th FORT, = (0 ( i), (2) , (3) ...0 (T+y) ) and 
oo = (oo(!), oo(2), oo(3)...oo(7' + y)), where a bracketed 
subscript enumerates the elements of or oo. We adopt 
a notation such that a sum of the form E$=oj where 
X and Y are the E-component vectors (Xi, X2 ■ ■ ■ Xp) 
and (Yi, Y% . . . Yp) respectively, denotes the set of sums 

Ex 1= o El=o Ex 3 =o • ■ ■ Ex' =o- Furthermore, we also 
assume that a state of the form \X) denotes the state 
\X\) |A 2 ) ■ ■ - \Xp). Note that due to an exchange sym- 
metry for atoms in the same group of atoms, it is suffi- 
cient to just consider the reduced density operators de- 
noted by Eqn l(2*T)l to deal with all possible pQ j 's. That 
is, we do not need to consider, say, tracing over the first 
and third atoms in the first FORT and the second one 
in the second FORT. This is so as the pQ j this yields is 
identical to that produced by tracing over the first two 
atoms in the first FORT and the first one in the second 
FORT. 

We now find [/OqJ 2 and then trace over the remaining 



2N — (T + V) atoms, producing 

(66 66 00,00,00,00 
EE E 
P=0 V=8 AfjA/ - ' ,M ,M ' 

f(Af,M')f(M,M') 
x (P\N)\N)(N'\{N'\P) 
x (V\M)\M)(M'\{M'\V)^ , (22) 

where, in analogy with P, V is a dummy variable given 

byP = (Pf) ...Vip ,V[ 2) -^v ) ) where V » ] denotes 
a vibrational number state in the x direction for the a th 
atom in the j th FORT. 

Using Eqn (J23, we now numerically determine 
Sl {PQj ) f° r arbitrary T and V particular values of N 
and r. Our results provide lower bounds for S{pq j ) 
as 5i(/9)/log 2 e < S(p) as can be verified by consider- 
ing a power series expansion for S(p). Hence, knowing 
SL(pQj) for all bipartite splits of IV'Cm) allows us to in- 
fer a lower bound for min(S a n) and hence one for Embe- 
We thus calculate all ^(pq-, ) for N = 2, 3,4 for a range 
of r values numerically using straightforward C++ code. 
These results are then used to place lower bounds on 
Embe ( | V'cm ) ) for 4-way, 6- way and 8-way entanglement 
which appear in Figs^Ja) and (b). 

As IV'Cm) is the sum of an infinite number of statevec- 
tors, to calculate Sl in practice, we truncate the sum 
over Af in the definition of \iPcm) at a finite value. This 
induces errors in our lower bounds for Embe (| "0cm)) for 
which upper bounds can be derived. For all data points 
in Figs 01 (a) and (b), the errors on our lower bounds for 
Embe(IV'cm)) have been calculated to be less than 10~ 3 
and hence are negligible. 

Two interesting features of Figs0](a) and (b) are that, 
firstly, for a given r value our lower bound on Embe 
decreases for increasing N. It is possible that we can un- 
derstand this behaviour by observing that for constant r 
we initially have a fixed entanglement resource, namely 
the entangled output of the NOPA. It is conceivable that 
the decrease under consideration results from this fixed 
resource being spread amongst a larger number of sub- 
systems as we increase N thus, perhaps, causing it to 
distribute less bipartite entanglement to any given bi- 
partite split of IV'Cm)- In turn, this may decrease the Sl 
of both halves of an arbitrary split, thus explaining the 
decrease in our lower bound for Embe for increasing N. 
The second interesting feature of Figs 01 (a) and (b) is 
that as we increase r, Embe increases as expected given 
that an increased r means that we have more centre-of- 
mass entanglement. 

V. DISCUSSION 

Throughout the paper, we have emphasised that IV'Cm) 
contains 2./V-way entanglement. However, for this entan- 
glement to be meaningful it must have observable effects. 
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FIG. 4: Plots of lower bounds on Embe (dimensionless) for 
|i/>cm) as a function of r (dimensionless) for a) N = 2 (4- 
way entanglement), b) (i) N = 3 (6- way entanglement) and 
(ii) TV = 4 (8-way entanglement). Note that in all figures we 
have linearly interpolated between points 0.1 units apart on 
horizontal axes. Numerical errors are less than 10 -3 for all 
data points. 



One feature of the system under consideration that makes 
its entanglement conducive to producing such effect is the 
fact that the atoms in the system are spatially separated 
and thus, in principle, individually accessible. Thus, for 
example, we could shine a sufficiently narrow laser beam 
on one of the atoms and, provided it did not propagate 
perpendicular to the x-axis, implement a local displace- 
ment on the vibrational state of the atom in the x direc- 
tion. Furthermore, accessing individual atoms is made 
easier by the fact that neighbouring atoms do not have 
be located at successive cavity-field nodes. Instead, they 
can occupy every second, third etc. node, thus increasing 
their spatial separation and making it easier to address 
them one at a time. Another advantageous consequence 
of the fact that each atom is individually accessible is 
that it permits us to perform measurements on the vi- 
brational states of single atoms, perhaps by employing 



a quantum-optical technique used to measure the posi- 
tion of individual trapped atoms by having them interact 
strongly with a low-photon number cavity mode |40| . 

In light of the considerations of the previous para- 
graph, some possible applications of the entanglement 
in IV'cm) are as follows: 

1. Violations of inequalities based on local realism. A 
number of such inequalities for an arbitrary num- 
ber of quantum systems have been formulated [3^| . 
Given the close connection between violations of 
these inequalities and entanglement, |i/'cm) is the 
sort of state we might expect to violate at least 
some 2iV-party inequalities based on local realism. 
However, as the Hilbert space for the vibrational 
motion each atom is infinite-dimensional and not 
two-dimensional as is the case for qubits, the vio- 
lations may require discretising or "binning" mea- 
surement results of a continuous variable such as 
quadrature phase amplitude. 

2. Solving quantum communication complexity prob- 
lems (distributed quantum computing). Quantum 
communication complexity problems |39| involve 
a number of parties attempting to evaluate some 
function / for a particular input string. Each party 
is given part of the input string and then uses 
shared prior entanglement, local classical compu- 
tation and public communication in attempting to 
evaluate /. In such a scenario, the prior entangle- 
ment can allow the evaluation to be performed in a 
superior manner to that attainable classically. As 
the entanglement in \ip)cM is such that every atom 
in the corresponding system is with every other one, 
it is a quantum resource seemingly well-suited to 
being of use to in solving quantum communication 
complexity problems better than can be done clas- 
sically. 

3. Continuous-variable quantum computation. 
Continuous-variable quantum computation |4l| 
involves quantum computing with infinite- 
dimensional quantum systems as opposed to the 
usual two-dimensional qubits. The most obvious 
way to perform this sort of computation with the 
system under consideration would be to, firstly, 
consider each atom in it as a qudit in the limit 
of d — ► oo. After this, we would then need to 
implement two-qudit gates by having different 
atoms interact with each other in a pairwise 
manner. One method by which this might be 
accomplished is by using a scheme |42T | employed 
in optical lattices to get two spatially-separated 
trapped neutral atoms of different species to 
interact with one another. This is done by varying 
the polarisations of the electromagnetic fields 
trapping the atoms which has the effect of varying 
the potentials that the atoms see in such a manner 
that they move towards each other. Once together, 
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the atoms interact via a dipole-dipole coupling. It 
is conceivable that this method could be applied 
to implement two-qudit gates in the system of 
interest. One complication, however, in utilising 
this scheme is that it necessitates that we modify 
our system by having the 2N atoms in it comprised 
of two different species, perhaps with the species 
of atom alternating as we move along each linear 
configuration. Nevertheless, whilst the system 
under consideration may not be the most natural 
one in which to do continuous-variable quantum 
computation, there is some possibility that the 
entanglement in it could be used to do this. 

A. Qualitative results 

The thinking underlying Definition 1 is the same as 
that which underlies the NPT sufficient condition for M- 
way entanglement. However, there are significant differ- 
ences between the two. Firstly, Definition 1 involves 
arbitrary dimensional subsystems, whereas the NPT suf- 
ficient condition deals only with qubits. Secondly, the 
NPT sufficient condition is a sufficient but not a nec- 
essary condition for M-way entanglement whereas the 
satisfaction of Definition 1 is both necessary and suffi- 
cient for pure states. Thirdly, the NPT sufficient condi- 
tion uses the partial transpose to determine the presence 
of M-way entanglement, whereas Definition 1 uses the 
mathematically simpler entity the trace of the square of 
a reduced density operator. Observe that Definition 
1 is narrower than the NPT sufficient condition in the 
sense that it only applies to pure states whilst the NPT 
sufficient condition is applicable to both pure and mixed 
states. 



to relate Embe to one or more quantum tasks or 
protocols such as distributed quantum computation 
with Embe telling us something valuable about how 
well we can perform these tasks. This is so as if we 
could do this, then it would increase Embe's util- 
ity. Unfortunately, however, this has not yet been 
accomplished. 

2. Can we tractably calculate Embe ? For an entangle- 
ment measure to be useful, it must be tractable and 
able to be calculated in practice. Unfortunately, 
Embe seems to be difficult to calculate, at least for 
the state considered. 

Although -Embe has the two above negative features we 
note that, firstly, further research may eliminate them 
and, secondly, we should consider them alongside the pos- 
itive features of Embe which are that it is a reasonable 
measure and that it helps us to understand the nature 
of the entanglement in |"0cm) and also the capabilities 
of quantum state exchange. Our results contribute to 
the understanding of multipartite entanglement involving 
massive particles and infinite-dimensional Hilbert spaces 
within a context that is not overly experimentally infea- 
sible. 

To conclude, we have shown that quantum state ex- 
change can be used to produce the state |^>cm) for two 
sets of trapped atoms in spatially separated FORTs. 
We have also show that IV'cm) is a 2 A- way entangled 
state and, in addition, have placed a lower bound on the 
amount of such entanglement that it possesses. Finally, 
we have discussed quantum information processing tasks 
that the 2A-way entanglement in |V'cm) could be used 
to help perform. 



B. Quantitative results 

A number of issues surround Embe , which we now dis- 
cuss. 

1. What does Embe tell us about what quantum re- 
source we have? Ideally, we would like to be able 
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